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Abstract 

I review the recent progress in measuring elliptic flow in heavy ion collisions. These measure- 
ments show clearly how hydrodynamics starts to develop as the system size is increased from 
peripheral to central collisions. During this transition, the momentum range described by hydro- 
dynamics increases as the system progresses from a kinetic to a hydrodynamic regime. Many of 
the systematic deviations from ideal hydrodynamics are reproduced effortlessly once the shear 
viscosity is included. In order to extract the shear viscosity from the data, kinetic theory can be 
used to determine which aspects of the elliptic flow reflect the details of the microscopic interac- 
tions, and which aspects reflect the underlying transport coefficients. I also review the identified 
hadron elliptic flow and the predictions of hydrodynamics for the LHC. 



1. Overview 

Perhaps the most important result from the Relativistic Heavy Ion Collider is the observation 
of strong elliptic flow (UEl- Elliptic flow is an asymmetry of particle production with respect 
to the reaction plane and has been measured as a function of transverse momentum, rapidity, 
and particle type. The interpretation of the observed flow which has been adopted by the heavy 
ion community is that the elliptic flow is the hydrodynamic response to the collision geometry. 
Implicit in this interpretation of the observed flow is that the time scale for momentum relaxation 
near the QCD phase transition is of order the quantum time scale 

h 

tr-—. (1) 

This estimate for the relaxation time is best expressed in terms of the shear viscosity to entropy 
ratio. For instance, in the viscous Bjorken model the energy density at time r evolves as 

de_e + P^4rj 
dr r 3rl' 

where e is the energy density, V is the pressure, and r] is the shear viscosity (3). Comparing the 
size of the viscous term to the ideal term, we conclude that hydrodynamics will provide a good 
description of the observed flow when 

« 1 . (3) 
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Using e + V - sT, and an estimate for the temperature and r , this criterion reads 



/iy/j\ /lftn\/ 300MeV 



) 
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From this estimate we see that hydrodynamics will begin to be a good approximation for rj/s< 0.3 
or so. To reiterate, implicit in the hydrodynamic interpretation of the flow results is a strong 
conclusion about the transport coefficients of QCD. 

Many of the systematic trends seen in the elliptic flow data do support the hydrodynamic 
interpretation of the observed flow. For example, as a function of centrality and transverse mo- 
mentum, the measured elliptic flow deviates from ideal hydrodynamics in a way characteristic of 
viscosity. These experimental trends have been fully clarified only recently and the experimental 
analysis is now rather sophisticated. These developments were reported on at the Quark Matter 
conference and are reviewed in Section [2j The systematics of the recent flow measurements give 
confidence in the overall picture of the hydrodynamic expansion. 

Although these experimental trends support the notion of a hydrodynamic response, there 
are several puzzling patterns in the elliptic flow data. For example, new data on the elliptic flow 
of the meson and the Q" baryon are reviewed in Section [2] The differences in the measured 
flow between mesons and baryons is generally explained with the coalescence model, which 
enjoys considerable phenomenological success. However, the coalescence model is theoretically 
unsatisfactory and is difficult to realize in a dynamical model. The seemingly simple coalescence 
trends seen in the elliptic flow data must be understood before the shear viscosity and other 
transport properties can be reliably extracted from the heavy ion data. 

In addition to experimental progress, there has been substantial theoretical progress in classi- 
fying the form of viscous corrections, both with viscous hydrodynamics and with kinetic theory. 
A brief summary of some of the developments discussed at Quark Matter 2009 is presented 
in Section [3] Most of these ideas presented in this summary are reviewed more completely in 
Ref.(2|, which was written shortly after the Quark Matter conference. Some sections from this 
longer review have been copied for this brief summary. 

2. Measurements 

One of the best ways to test the hydrodynamic interpretation is to systematically observe 
how the response changes from small systems to large systems. Experimentally this can be 
accomplished by colliding small nuclei such as CuCu and selecting peripheral collisions. Unfor- 
tunately measuring elliptic flow in these smaller systems is rather difficult, and reliable, precise 
results have become available only recently. 

The difficulty in measuring flow in small systems stems from fluctuations. Especially in 
peripheral AuAu and CuCu collisions, there are fluctuations in the initial eccentricity of the 
participants. Thus rather than using the continuum approximation to categorize the geometry, it 
is better to implement a Monte-Carlo Glauber calculation and estimate the eccentricity using the 
"participant plane eccentricity". This event by event eccentricity is denoted e PP in the literature. 
Clearly the experimental goal is to extract the response coefficient C relating the elliptic flow to 
the eccentricity on an event by event basis 



V2 = Cepp . 



(5) 
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Figure 1: (a) The elliptic flow measured by different methods |6|. The methods measure different quantities and should 
agree only when fluctuations are neglected, (b) The average elliptic flow <V2> measured by the different methods. To 
determine <V2> for each method, non-flow was estimated using measured pp data and fluctuations were estimated with a 
geometric model. 



If the flow methods measured (v2>, then we could simply divide the measured flow to de- 
termine the response coefficient, C = (v2> / (epp). The PHOBOS collaboration deciphered the 
confusing CuCu data by recognizing the need for e PP and following this procedure[4|. However, 
it was generally realized (see in particular. Ref.Q) that the elliptic flow methods do not measure 

precisely (v2>. Some methods (such as two particle correlations V2 {2}) are sensitive to J{^j), 
while other methods (such as the event plane method V2 {EP}) measure something closer to (V2). 
What precisely the event plane method measures depends on the reaction plane resolution in 
a known wayO. So just dividing the measured flow by the average participant eccentricity is 
not entirely correct. The appropriate quantity to divide by depends on the method O El El - In 
a Gaussian approximation for the eccentricity fluctuations this can be worked out analytically. 

For instance, the two particle correlation elliptic flow V2 {2} (which measures -My^j), should be 

divided by ^J^^p~p^) . With a complete understanding of what each method measures, Ref.|6j was 
able to make a simple model for the fluctuations and non-flow and show that (v2> measured by the 
different methods are compatible to high precision. This is illustrated in Fig.[T] The stunning pre- 
cision of the recent flow measurements will be very useful in understanding how hydrodynamics 
begins to work in peripheral collisions. This work should be extended to the CuCu system where 
non-Gaussian fluctuations are stronger and ultimately corroborate the PHOBOS analysis EJ [9l- 
This is a worthwhile goal because it will clarify fully the transition into the hydrodynamic regime 

EQlCLO. 

The scaled elliptic flow V2/6 measures the response of the medium to the initial geometry and 
can be used to determine the regime of validity of hydrodynamics. Additional information about 
the shear viscosity can be gleaned from the transverse momentum dependence of the observed 
elliptic flow. Fig. [2] shows V2(pr)/£ as a function of centrality, 0-5% being the most central and 
60-70% being the most peripheral. Examining this figure we see a gradual transition from a weak 
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Figure 2: Elliptic flow V2(pr) as measured by the STAR collaboration (T^E) for different centralities. The measured 
elliptic flow has been divided by the eccentricity. The curves are ideal hydrodynamic calculations based on Refs. ll4l|75l 
rather than the viscous hydrodynamics discussed in much of this summary. 

to a strong dynamic response with growing system size. The interpretation adopted here is that 
this change is a consequence of a system transitioning from a kinetic to a hydrodynamic regime. 

There are several theoretical curves based upon calculations of ideal hydrodynamics fT5l [T6l 
which for p T < 1 GeV approximately reproduce the observed elliptic flow in the most central 
collisions. Since ideal hydrodynamics is scale invariant (for a scale invariant equation of state) 
the expectation is that the response V2/6 of this theory should be independent of system size or 
centrality. This reasoning is borne out by the more elaborate hydrodynamic calculations shown 
in the figure. On the other hand, the data show a gradual transition as a function of increasing 
centrality, rising towards the ideal hydrodynamic calculations in a systematic way. These trends 
are captured by models with a finite mean free path lfTTl . 

The data show other trends as a function of centrality. In more central collisions the linearly 
rising trend, which resembles the ideal hydrodynamic calculations, extends to larger and larger 
transverse momentum. Viscous corrections to ideal hydrodynamics grow as 

where L is a characteristic length scale. Thus these viscous corrections restrict the applicable 
momentum range in hydrodynamics fTjfl . In more central collisions, where £ mfp /L is smaller, 
the transverse momentum range described by hydrodynamics extends to increasingly large p?. 
These qualitative trends are reproduced by the more involved viscous calculations Q. 

While many of the trends seen in Fig. [T] and Fig. [2] are reproduced and understood with 
viscous hydrodynamics, there are additional trends in the elliptic flow data which are only par- 
tially understood. For instance Fig. [3Ja) shows the elliptic flow of identified particles 7r, K, p. 
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Figure 3: (a) The elliptic flow as a function of transverse momentum for identified particles as measured by the PHENIX 
collaboration fl9l . (b) The elliptic flow of identified hadrons rescaled according to the quark coalescence model and 

plotted as a function of KEj = Jp\ + M 2 - M. 



At low momenta the separation amongst the different particle species is well reproduced by hy- 
drodynamics. However as the momentum is increased the proton elliptic flow equals, and then 
exceeds, the pion elliptic flow. These systematic trends are seen in all collision systems and 
centralities. The prevailing explanation is that constituent quarks coalesce to form hadrons at the 
phase boundary. This ansatz is supported by the observation that if the hadron momentum and V2 
is divided by the valence quark content all of the V2 of the different hadron species lie along a sin- 
gle curve. This is illustrated in Fig.^b) which is plotted as a function of KE T = ^Jp* + M 2 - M 
rather than transverse momentum to capture the hydrodynamic behavior at small momentum. Al- 
though constituent quark scaling works rather well, the theoretical support for quark coalescence 
is small since it is difficult to realize a coalescence mechanism in a dynamical simulation. It 
nevertheless remains to find an alternative picture for the observed different flows of mesons and 
baryons. At the Quark Matter conference the elliptic flow of identified particles was measured 
accurately out to rather large transverse momenta. At sufficiently large momenta the data deviate 
from a universal coalescence curve providing new insight into the hadronization dynamics in this 
region. 

To conclude this section, we turn to Fig. [4] which compares the elliptic flow protons and pions 
to the flow of the multi-strange hadrons Q" and 0. The important point is that the Q" is nearly 
twice as heavy as the proton and more importantly, does not have a strong resonant interaction 
analogous to the A. For these reasons the hadronic relaxation time of the £l~ is expected to be 
much longer than the duration of the heavy ion event. Nevertheless the £l~ shows nearly the 
same elliptic flow as the protons. This provides fairly convincing evidence that the majority of 
the elliptic flow develops during a deconfined phase which hadronizes to produce a flowing Q" 
baryon. 
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Figure 4: A comparison of the elliptic flow of pions and protons to the elliptic flow of the multi- strange and Q f20l . 



3. Modeling Elliptic Flow with Hydrodynamics and Kinetic Theory 

Many of the trends seen in the data (with the possible exception of quark coalescence) are 
reproduced by hydrodynamics and kinetic theory. Generally the kinetic theory estimates for the 
shear viscosity are consistent with the estimates from viscous hydrodynamics. Specifically unless 
rj/s<03 it is impossible to reproduce the observed elliptic flow. Given that the transport time 
scales extracted from the heavy ion data are close to the quantum time scale given in Eq. ([I]), 
it is clear that the microscopic details of kinetic models cannot be trusted. Nevertheless these 
shortcomings of the microscopic theory are unimportant in the hydrodynamic regime. In the 
hydrodynamic regime the only properties that determine the evolution of the system are the 
equation of state, Pie), and the shear viscosity and bulk viscosities, r\(e) and f (e). In the sense 
that kinetic theory provides a reasonable guess as to how the surface to volume ratio influences 
the forward evolution, these models can be used to estimate the shear viscosity, and the estimate 
may be more reliable than the hydrodynamic models. More importantly, by comparing the results 
of different microscopic models one can determine which features of the heavy ion data are 
universal (i.e. only depend on 77(e), V(e) and g(e) ). 

There were a number of promising efforts to reproduce hydrodynamic results from kinetic 
theory reported at the conference. First there was an effort to reproduce hydrodynamic shocks 
with kinetic theory by the Frankfurt group. Fig. [5ja) shows a kinetic theory simulation of the 
shock tube problem. Clearly the BAMPS code is capable of reproducing the correct hydro- 
dynamic limit in detail. The deviations of the BAMPS code from ideal hydro are beautifully 
reproduced by viscous hydro. This gives a great deal of confidence in the BAMPS code and in 
the viscous hydrodynamic code vSHASTA. Additional results from the BAMPS simulation of 
elliptic flow were presented in the poster session. 

A similar approach to the hydrodynamic limit was reported by Huovinnen and Molnar |2T1L 
and Gombeaud |11]. In particular Fig. [5jb) shows a simulation of the MPC code which also 
makes a direct comparison with viscous hydrodynamics for the Bjorken expansion of an ideal 
massless gas with constant cross section. As the Knudsen number K = cr/nR 2 dN/dy is in- 
creased, the simulation approaches Israel-Stewart hydro and ultimately the Navier-Stokes limit. 
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Figure 5: (a) The pressure relative to the initial pressure for shock tube initial conditions. The three curves are for the 
BAMPS parton cascade model, viscous hydrodynamics (vSHASTA), and ideal hydrodynamics [22). (b) The ratio of 
longitudinal and transverse pressures in kinetic theory, Israel Stewart hydrodynamics, and the Navier-Stokes equations, 
for a Bjorken expansion of a massless ideal gas with constant cross section. The simulations are compared as a function 
of Knudsen number, K = cr/nR 2 dN/dy (21). 

Given that the kinetic code and the viscous hydrodynamic simulations agree reasonably, the MPC 
code can be used to reliably extract the shear viscosity from the heavy ion data. 

When viscous hydrodynamics is extended to a second order there are additional relaxation 
times (e.g. r n ) beyond the shear viscosity, rj(e). Just as transport models should be approximately 
independent of the details of the microscopic interactions, results from viscous hydrodynamics 
should be approximately independent of the precise way in which the second order terms are im- 
plemented. This is indeed the case (23H2UG3. Additional results from viscous hydrodynamics 
will be discussed more completely by P. Romatschke in this volume [26]. 

4. Outlook 

Clearly there is a strong convergence between kinetic and hydrodynamic simulations of heavy 
ion reactions. These simulations reproduce many trends observed in increasingly precise mea- 
surements of elliptic flow. This convergence strongly suggests that the hydrodynamic interpre- 
tation of the observed flow is correct. One of the striking tests of hydrodynamic predictions is 
the saturation of elliptic flow at high energy. From RHIC to the LHC, hydrodynamics predicts 
an increase in the flow which is significantly less than a naive extrapolation from lower energies. 
This is illustrated in Fig. [6] and will be one of the first tests of the hydrodynamic paradigm at the 
LHC. 
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Figure 6: Figure presented by W. Busza providing an estimate for the elliptic flow at the LHC by extrapolating trend 
lines from lower energy to higher energy (solid black lines) [27 1. The dashed red lines shows an estimate based on ideal 
hydrodynamics (28). 
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